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Abstract 

Determining the worst-case uncertainty added by a quantum circuit is shown to be compu- 
tationally intractable. This is the problem of detecting when a quantum channel implemented 
as a circuit is close to a linear isometry and it is shown to be complete for the complexity 
class QMA of verifiable quantum computation. The main idea is to relate the problem of de- 
tecting when a channel is close to an isometry to the problem of determining how mixed the 
output of the channel can be when the input is a pure state. 

1 Introduction 

A linear isometry U: "K — > % is a linear map that preserves the inner product of any two ele- 
ments, or equivalently satisfies U*U = t^. These transformations are fundamental in quantum 
computation: they are exactly the maps that may be realized using unitary quantum circuits with 
access to ancillary qubits in a known pure state — the standard model of quantum computation. 
It is an important problem to determine when a computation in a non-unitary model, such as 
measurement based quantum computing or computation in the presence of noise, approximately 
implements some operation in the unitary circuit model. In this paper it is shown that this problem 
is QMA-complete when the input computation is modelled as a quantum circuit consisting of the 
usual unitary gates, plus the ability to discard qubits as well as introduce ancillary qubits. The 
circuit model is not essential: the hardness result also applies to any model that can efficiently 
simulate and be simulated by the mixed-state circuit model. 

The complexity class QMA is the quantum analogue of N P: the class corresponding to classically 
verifiable computation. This concept was first considered in fill], first formally defined in [a], and 
first studied in jl^l . QMA is the class of all problems that can be verified with bounded error 
by a polynomial-time quantum verifier with access to a quantum proof. This proof is given by a 
quantum state on a polynomial number of qubits and may depend on the input. 

The class QMA has complete (promise) problems: problems in QMA that are computationally 
at least as hard as any other problem in the class. This implies that an efficient algorithm for any 
of these complete problems can be used to find an efficient algorithm for any problem in QMA. 
The simplest of these complete problems is the 2-local Hamiltonian problem, which is informally 
the quantum version of the circuit satisfiability problem for unitary circuits with gates of constant 
size. A formal description of this problem, as well as a proof that the 5-local Hamiltonian problem 
is QMA-complete can be found in Jioll . The improvement of this result to the 2-local case is due 
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to Kempe, Kitaev, and Regev 
local consistency H3] (see also 




Several other complete problems for QMA are known, such as 
some problems related to the minimum output entropy | 



testing whether unitary circuits are close to the identity |@] (see also |@]), and finding the ground 
states of some physical systems (lallZD- In the present paper we add a new complete problem to 
this list: the problem of determining if a quantum circuit implements an operation that is close 
to an isometry. As discussed in Section [3J this is equivalent to determining if the channel always 
maps pure states to states that are approximately pure. 

The remainder of the paper is organized as follows. Section |2] introduces notation and back- 
ground. Section|3]introduces the notion of approximate isometries and makes formal the problem 
of detecting when a channel is an approximate isometry. The QMA-hardness of this problem is 
proved in Section|4]and proof of the containment in QMA, the most technical portion of the result, 
appears in Section [5] 



2 Preliminaries 

In this section the notation and background that is used throughout the paper are presented. Much 
of the notation used here is standard and this is in no way a complete introduction to quantum 
information. See lfl5ll for a more detailed treatment of these topics. 

All Hilbert spaces considered in this paper are assumed to be finite-dimensional and are 

denoted by scripted capital letters "K, X, The pure states are the unit vectors in these spaces. 

The set of density matrices or mixed states on "K is given by D(!H), and the set of all quantum 
channels mapping D(!K) to D(3C) is T["K,X). The quantum channels are exactly the completely 
positive and trace preserving linear maps. The identity channel in T(IK,IK) is denoted while 
1 -ji is the identity on 'K. 

Given a quantum channel <D G T(5f, X) we make use of two representations. The first of these 
is the Choi representation Jj], which provides a unique representation of a channel O € T(!K, X) 
as a linear operator on X "K. This representation is given by C(O) = (O <g> I^)(I4> + )(<|) + |), where 
|c{) + ) = Y_i fti) / is a maximally entangled state in "K <S> "K. 

The second representation that we use is the representation of a completely positive map O by 
a set of Kraus operators: matrices At such that ® (X) = Y-i AiXA? . This representation is also due 
to Choi Q]. If in addition the map O is trace preserving, then the operators At satisfy the property 
Y_i A* Ai = 1 . The number of Kraus operators in a minimal Kraus decomposition is given by the 
rank of the Choi matrix C(O). 

In order to measure how close a state is to being pure we use the operator norm HX^, which 
for a linear operator X is the largest singular value of X. When X is normal, this is simply the 
largest eigenvalue (in absolute value) of X. Dual to the operator norm is the trace norm, which for 
a linear operator X is given by ||X|| tr = tr \/X*X. This is exactly the sum of the singular values of 
X. When X is a quantum state, this simplifies to the sum of absolute values of the eigenvalues of 
X, so that || p || tr = 1 for all density matrices p. 

One final quantity that we use is the fidelity, which for two density matrices p, cr is given by 
F(p,ff) = tr yj y/pa y/p. While it is not obvious from this definition, the fidelity is symmetric 
in the two arguments. When one of the arguments is a pure state, the fidelity simplifies to 
F(p, |i[>)(i[>|) = yf (ip|p|i|;). An important relationship between the trace norm and the fidelity is 

2-2F(p,W;><iH) 2 < llp-I^MIItr (1) 
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that we use to relate different notions of the purity of a quantum state. This inequality can be 
found in |fl5l. Chapter 9]. 

We require one final piece of background. In order for a quantum channel to be given as 
input to a computational problem we need a representation of the channel. Using either the Choi 
matrix or Kraus operators produces a representation that, in the case of channels implementing 
efficient quantum algorithms, is exponentially larger than the size of a circuit representation. These 
channels have circuit representations that are logarithmic in the Hilbert space dimension. For this 
reason, we use a circuit representation of quantum channels. Such a representation is provided by 
the mixed-state circuit model of Aharonov et al. which is simply the usual model of unitary 
quantum circuits with two additional gates. These gates are the gate that introduces ancillary 
qubits in the |0) state and the gate the traces out (i.e. discards) a qubit. This circuit model can be 
used to represent any quantum channel, which makes it ideal for the problem that we consider. 



3 Isometries and rank non-increasing channels 

One important property of the linear isometries is that they do not increase rank. This is essential 
to the QMA protocol in Section |5l which is able to detect exactly those channels that are rank- 
increasing. More formally, a channel O is rank non-increasing if for all states p the output of (D 
satisfies rank(p) ^ rank(0(p)). Unfortunately, this property does not characterize the isometries. 
Consider the channel O ( p) = |0) (0| that discards the input state and returns a fixed pure state. This 
channel is not an isometry but it is also rank non-increasing. 

This property can be used to characterize the isometries if we make a small adjustment. The 
channels that are rank non-increasing when adjoined to an auxiliary space of arbitrary dimension 
are exactly the isometries. We call a channel O G T(IK, X) completely rank non-increasing if for 
any 1 the channel O <S> l-j is rank non-increasing, i.e. if rank[(0 ® Ij-)(p)] ^5 rank(p) for all p. 
The channel <D(p) = |0)(0| is not completely rank non-increasing: consider applying it to half of a 
maximally entangled state (O <g> I:n)(]c|> + )(ct> + |) = |0)(0| <g> Ijc/ dim "K. As in the case of complete 
positivity, we need only to verify this property on an auxiliary space of the same dimension as the 
input space. It is also easy to see that this property characterizes the linear isometries. 

Proposition 1. The following are equivalent for a channel O G T(!K,3C): 

1. ®(p) = UpU*/or some linear isometry Mfrom "K to X, 

2. O is completely rank non-increasing, 

3. O ® Lk is rank non-increasing. 

Proof. The first two implications are immediate. To prove that © =4> CD, let O ® 1^ be rank 
non-increasing. This implies that rank(C(0)) = 1. Recalling that the number of Kraus operators 
in a minimal decomposition is rank(C(0)), it follows that O can be expressed as O(p) = ApA*. 
The condition that O is trace preserving implies that the operator A satisfies A* A = tjt- □ □ 

This characterization guides the remainder of the paper. Detecting when the channel <D <g> 1^ 
increases rank provides an operational method to determine when a channel is an isometry. 
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3.1 Approximately pure states 

In order to show that non-isometry detection is QM A-complete we need to consider an approximate 
version of the problem. This is because a protocol for a QMA language is permitted to fail with 
small probability. The definition of approximate isometries used here is closely related to the 
notion of approximately pure states. Several equivalent notions of the purity of a density matrix 
are considered in this section. 

Perhaps the most well-known notion of how close a mixed state p is to being pure is the purity 
of p, given by tr(p 2 ). A similar measure is given by || pH^, the largest eigenvalue of p. It is not 
hard to see that these quantities are related. If p = 2Ii Ai|Ai)(A-jJ is the spectral decomposition 
of p, with the eigenvalues in decreasing order, then tr p 2 = ^ i A? ^ A 2 = || pH^. In the other 
direction, since the purity is convex, it is maximized for 1/Ai eigenvalues each of value \\, i.e. 
tr p 2 = Y-\ ^? ^ A 2 /Ai = || P lloo- Taken together, these two inequalities show that 

llpll^tr^KHplL. (2) 

These quantities are also related to the more familiar trace distance on quantum states. 

Proposition 2. Let p e D(!K) and let e > 0. There exists a pure state € "Ksuch that || p — |4>)(4>||| tr ^ 
e ifandonly zf HpH^ ^ 1 — e/2. 

Proof. Let p have spectral decomposition given by p = Y.\ Ai|Ai)(At|, with A} ^ A2 ^ . . . ^ A^. If 

II P lloo = A i > 1- e/2, then 

||p-|A 1 )(A 1 ||| tr = (l-A 1 ) + ^ d =2 A i = 2(l-A 1 ) <2(e/2) = e. 

On the other hand, if |i^) G "K is a state such that || p — |ip) (i|)| || tr e, then by Equation ((TJ) 

e > || p - W»> (M || tr > 2 - 2 F(p, \^)(^\) 2 = 2- 2(^|pW>) = 2 - 2^. A t | {M\) | 2 . (3) 

The final quantity is a convex combination of the A^, with weights determined by the state |i[)). 
This is maximized when |i|>) = |Ai), since Ai is the largest eigenvalue of p. Combining this with 
Equation (0 we have e ^ 2 — 2Ai = 2 — 2 || p H^, which implies that || p ^ 1 — e/2. □ □ 

Given these notions of purity, we will call a state z-pure if || p ^ 1 — e. By the previous results 
the purity of such a state satisfies tr(p 2 ) ^ (1 — e) 2 ^ 1 — 2e, and there is a pure state |ip) such that 
|| p — |ip) (4>| || tr ^ 2e. For the results of this paper, any of these three measures suffices, as they are 
equivalent up to polynomial factors in e. 

3.2 Approximate isometries 

The focus of this paper is to show that detecting when a channel is far from an isometry is 
computationally difficult. To do this we need to define the class of channels that are the approximate 
isometries. Isometries always map pure states to pure states, even in the presence of a reference 
system. Proposition [TJ shows that this condition characterizes the isometries. Weakening this 
requirement, we call a channel an e-isometry if it maps pure states (over the input space and a 
reference system) to states that are e-pure, for some e > 0. 
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More formally a channel cp g T(IK,3C) is an e-isometry if for any pure state |i|>) G IK £g> K the 
output of O Ijc satisfies || (O Ijc)(I 1 ^)(' 1 H) II oo ^ 1 — £ / i- e - when applied to part of any pure state 
the output state is close to pure. This implies that O Iji does not reduce the operator norm of 
any input by more than a factor of 1 — e. We use this to define the computational problem that is 
the main focus of the paper. 

Problem 3 (Non-isometry). For ^ £ < 1/2 and a channel <D G T(IK, DC), given as a mixed-state 
quantum circuit, the promise problem is to decide between: 

Yes: There exists a pure state G IK such that || (O FhHKWW^I) < £, 

No: For all pure states |it>) G IK, ||(<D I;k)(N>)M>I)||oo > 1 ~ £ - 

When the value of e is significant, we will refer to this problem as Non-isometry £ . 

Using the equivalence results of Equation (|2j) and Proposition |2j this problem may be equiv- 
alently defined in terms of either the purity or the trace distance to the closest pure state, up to 
a small increase in e. The case of the minimum output purity of a channel has been studied in a 
different context by Zanardi and Lidar l2oll , though they focus on finding the minimum purity of 
a channel over a subspace of the inputs. The problem we consider here is equivalent to evaluating 
the channel purity of <D 1^ over the whole input space. 

The difficulty of the Non-isometry problem does not change if the dimension of the ancillary 
system is permitted to be larger than the size of the input system, so long as the number of qubits 
needed to represent the ancillary system is polynomial in the number of input qubits. 

The notion of approximate isometry that we consider here is not equivalent to the channel being 
completely rank non-increasing on average. This property is modelled by the distance between 
the Choi matrix of a channel and a pure state. While it is true that the Choi matrix is pure if and 
only if the channel is an isometry, it is close to pure in the trace distance when the channel is close 
to an isometry on average. In this paper we consider the worst-case, i.e. we consider a channel to be 
close to an isometry if and only if the output of <P £g> I •x is close to pure for any pure state input. A 
simplification of the protocol presented in Section [5] yields a polynomial-time quantum algorithm 
for the problem of determining how close the Choi matrix of a channel is to a pure state. This is 
because C(O) can be generated efficiently, and given two copies the swap test can be used to test 
the purity of a quantum state as shown in jjj. 

4 QMA hardness 

In order to prove the hardness of Non-isometry we modify an arbitrary QMA protocol to obtain a 
circuit that can output a mixed state exactly when the verifier would have accepted in the original 
protocol. This yields a circuit that is far from an isometry if and only if there is a witness that causes 
the verifier in the original protocol to accept. Deciding whether or not there is such a witness is 
QMA-hard, by the definition of the complexity class. More formally, a language L is in QMA if 
there is a quantum polynomial-time verifier V such that 

1. if x G L, then there exists a witness p such that Pr[V accepts p] ^ 1 — e, 

2. if x i L, then for any state p, Pr[V accepts p] < e, 
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Figure 1: Verifier's circuit in a QMA protocol. The verifier accepts the witness state if and only 
if the measurement in the computational basis results in the |1) state. 
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Figure 2: Constructed instance of Non-isometry. The output state is mixed by the completely 
depolarizing channel O only if the state is a valid witness to the original QMA protocol. 

The exact value of the error parameter e is not significant: any e < 1/2 that is at least an inverse 
polynomial in the input size suffices jlol,E3l . 

Let L be an arbitrary language in QMA, and let x be an arbitrary input string. The goal is to 
embed the QMA-hard problem of deciding if x G L into the problem of testing whether a mixed- 
state quantum circuit is close to an isometry. Let V be the isometry representing the algorithm of 
the verifier in a QMA protocol for L on input x. We may "hard-code" the input string x into V 
because the circuit needs only to be efficiently generated from x. The algorithm implemented by 
the verifier is shown in Figure [U The verifier first receives a witness state |i[>), applies the isometry 
V, and then makes a measurement on one of the qubits, the result of which determines whether 
or not the verifier accepts. Any qubits not measured are traced out. 

For concreteness, let V act on the input spaces W and A, which hold the witness state and the 
|0) state of the ancilla respectively. Let M be the space corresponding to the measured output qubit 
in the protocol and let 3 represent the 'garbage' qubits that are traced out at the end of the protocol. 
The probability that verifier accepts the witness state € W is 

Pr[V accepts M>>] = (l|tr s [V(W>)(il>| ® |0>(0|)V*] |1). (4) 

Deciding if there is some such that this expectation is close to one is complete for QMA. 

From Figure[TJit is simple to construct a circuit that produces highly mixed output exactly when 
there exists such a The idea is add a controlled application of the completely depolarizing 
channel D. on the space 9/ instead of tracing it out. The resulting circuit is shown in Figure |2l In 
the case that the verifier accepts with negligible probability for every input state |i[>), then both 
the measurement and the controlled depolarizing channel have little effect, leaving the state of the 
system close to a pure state. If, on the other hand, there is a state on which the verifier accepts 
with high probability, then on this input the circuit in Figure |2] produces a highly mixed state. 
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Formalizing this notion proves that Non-isometry is QMA-hard. 

Theorem 4. Let e > be a constant, and let p be the maximum acceptance probability of the protocol V. 
Let G> e T(W,M S) fee f/ie cz'rcuz'f m Figure^ Then if dim 3? = dim W 

p<^ £ => min||(0®I K )(|^)(i|;|)|| 00 > 1 - e, 

p5sl-e min||(0® I K )(|ip)(iH)|| 00 < e. 

Proof. Notice that we may assume that the output dimension of O is dimM <2> 9 = 2d > 2/e by 
padding the circuit for V with log 1/e unused ancillary qubits, if necessary. 

As a first step, we evaluate the output state of the channel <D <g> 1^. Applied to a pure state 
|t|>) £ W®X this channel first adds the ancillary |0) qubits in the space A and then applies the 
isometry V from the QMA protocol. This results in the pure state |4>) = (V® ® |0)). We 

may decompose this state in terms of the qubit in the space M, obtaining for some ^ p ^ 1 

I4>) = 71^10) ® |4) ) + VpII) 

The value of p is exactly the probability that the measurement result is |1), i.e. the probability that 
the verifier will accept the input state tr^ |ol>)(i|)| in the original protocol. Using this, the state after 
the measurement and the controlled depolarizing channel on 9 is 

(1 - p)|0)<0| 1+oX+ol + (p/d)|l)(l| ® l s ® p, (5) 

where p is the residual state on 01 after this channel has been applied (p = trg I4>i}(4 5 iL but this will 
not be important). Evaluating the largest eigenvalue of this state we find that 

||(0®I 3 ,)(|4))(4)|)|| 00 =max{l-p,^||p|| 00 }. (6) 

We analyze the maximum in Equation ([6]) in two cases. The first of these cases is when there is 
no input the verifier accepts with probability larger than e. In this case the output of the channel 

(g> Igj is given by Equation (O where p ^ e. Here Equation |6]) shows that the output has an 
eigenvalue of magnitude at least min^ || (O £g> I^)(||x)((x|) || ^ > 1 — p > 1 — e. 

The second case is when there exists a state that verifier to accepts with probability at least 

1 — e. In this case we take the input state to O ® I# to be \y) = <g> |0), i.e. we set the reference 
system to be any pure state that is not entangled with the rest of the input. The output is given by 
Equation (O withp > 1 — e and p = |0)(0|. Equation ® yields 

miri||(0(8)I K )(|p)(M.|)|| 00 < ||(0®I 3l )(|Y)(y|)|| 00 =max{l-p,^ HpII^} < max je,^j = e, 

as we have taken 1/ d < e (by adding O (log 1/e ) unused ancillary qubits if necessary) . □ □ 

This theorem shows that determining how far the output O <8> Id? is from a pure state is as 
computationally difficult as determining whether or not the verifier can be made to accept with 
high probability in a QMA protocol. Since the construction of the circuit shown in Figure |2] can be 
performed efficiently, this implies the hardness of this problem. 

Corollary 5. For any constant < e < 1/2, Non-isometry is QMA-hard. 

Using the equivalences between notions of purity in of Section l3Tl this also implies that evaluating 
the purity of a quantum channel, as defined by Zanardi and Lidar ll^oh is QMA-hard. 
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5 QMA protocol 



In order to show that Non-isometry is QMA-complete, it remains only to construct a QMA protocol 
for the problem. The key idea behind this protocol is that when two copies of a channel <P are 
applied in parallel to the input state <g) the output lies in the antisymmetric subspace if and 
only if <D(|i|))(ip|) is a mixed state. This provides a probabilistic test that can detect when a channel 
is far from an isometry. 

Unfortunately, in a QMA protocol the verifier cannot assume the witness is given by two non- 
entangled pure states. It suffices, however, for the verifier to require that the input state lies in the 
symmetric subspace of the input space {% ® ft)® 2 . To show that the channel is not an isometry 
in QMA, the prover can provide a symmetric state that a parallel application of the channel maps 
into the antisymmetric space of the output space [X (g) ft)® 2 . 

The verifier in such a protocol needs a test to determine when a state is symmetric or anti- 
symmetric. Such a test is provided by the swap test, which was introduced in the context of 
communication complexity in U, though we make use of it to test purity using an idea from 

The swap test can be characterized as the projection onto the symmetric and antisymmetric 
subspaces of a bipartite space. If W is the swap operation on a space 5C ® !K, then the symmetric 
measurement outcome of the swap test corresponds to the the projector {tjt®3t + W)/2, and the 
projector (Ijcg^ — W)/2 corresponds to the antisymmetric outcome. 

The main idea behind the protocol for Non-isometry is that the swap test can be used to measure 
the purity of a state. As observed in [5j], when applied two to copies of a state p = Y.i ^il4>i)(4>il 
the swap test returns the antisymmetric outcome with probability 

\ tr((l - W)(p ® P)) = \ - = \ ~ \ MP 2 )- (7) 

This implies that the swap test on two copies of a state can be used to test purity and, by extension, 
when a channel is far from an isometry. 

A straightforward protocol for Non-isometry on a channel <D is then to receive a witness state 
W) ® W), a pply the channel to obtain [(O ® I)(|i[>)(il>|)]® 2 , and finally apply the swap test. The 
result is the antisymmetric outcome with high probability only when the state (<D ® I)(|i|>)(i|;|) is 
highly mixed. Such a protocol detects the channels that are far from isometries. 

Unfortunately, the verifier in a QMA protocol cannot assume that the witness state is of the form 
|i|>) <g> The verifier can check that he has received some state in the symmetric subspace and 
then use the fact that this subspace is closed under the parallel application of a rank non-increasing 
channel. The verifier in the following protocol uses the swap test to both check the symmetry of 
the input state and the antisymmetry of the output state. 

Protocol 6 (Non-isometry). On an input channel O £ T(5C, %): 

1. Receive a witness state p G D((5£ <g> ft)® 2 ), where ft is a reference space such that dimft = 
dim "K. Apply the swap test to p, rejecting if the outcome is antisymmetric. 

2. Use the channel O to obtain cr = (<D ® I^)® 2 (p). 

3. Apply the swap test to a, accepting if the outcome is symmetric and rejecting otherwise. 

A diagram of this protocol can be found in Figure |3l The correctness of this protocol is argued in 
the following theorem. 
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Figure 3: QMA protocol for Non-isometry. The verifier accepts only if the first swap test results 
the symmetric outcome and the second swap test results in an antisymmetric outcome. 

Theorem 7. Let O e T{!H,X), and let p(p) be the probability that the verifier described in Protocol® 
accepts the input state p e D((!K <S> H)® 2 ), then 

1. If mm^ || (O tg) H^, ^ e, fen fere exzsfs a witness p swc/z fef p(p) ^ (1 — e)/2. 

2. P/min^ ||(® <g> I^HM^^iM) ^ 1 — t, then for any witness p,y{p) ^ 9e. 

Proof. For the sake of brevity let <J> = O ® 1^ throughout. To prove the first assertion, let |ip) be 
a pure state in <8> 3? for which || ( I , (N>)( l l , l)||oo ^ £ / ar *d let the witness state p = I4>)(i[>| ® \ty)($>\. 
This state is invariant under the swap operation and so the swap test in Step H] passes and does 
not change the state. Step|2]results in the state [®(|ip)(iH)]® 2 . Using Equations 10 and 0, the final 
swap test returns the antisymmetric outcome with probability 

i - |tr [<f>(W>>U>|) 2 ] >\-\ ||*(H>>M>I)L > ^p, 

and so the verifier accepts p with probability approaching one-half for small e. 

To show the second assertion, we take 4> is an e-isometry and analyze the probability that 
the verifier can be made to accept. We may assume that the witness state lies in the symmetric 
subspace of ["K ® IR)® 2 , as the verifier either rejects in Step [T] or projects the witness onto this 
subspace. To complete the proof, we show that (<J>)® 2 leaves p approximately symmetric. 

To do this, we approximate O by an operator that preserves the symmetry of input states. Let 
(Ji) : 1 ^ i ^ dim IK} be an orthonormal basis for the spaces "K, IR (this is possible because they 
have the same dimension). The states {|tj) : 1 ^ i, j ^ dim Ji} are an orthonormal basis for 'K <g> 3?. 
Since <J> approximately preserves rank, there are states |i|h) G % such that 

|| (CD ® IaH|ij)(ij|) " NH>W>il ® UXJHItr < £ (8) 

for all i and j. We define a linear operator A: "K — > % by the equation A|i) = Cil4>i), where the 
Ci £ C with \c\\ = 1. The introduction of the phases Ct is necessary because Equation $8$ only 
defines the states |iJh) up to a phase. Note that the operator A is not necessarily unitary as we may 
not assume that the states are orthogonal. The next step is to show that, for some choice of 
the phases Ct, conjugation by A approximates the channel <D in the trace norm. This is the most 
technical portion of the proof. 
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Consider the output of O on the entangled state + )/ V2 in "K ® % given by 

p = i Y. *daa)(bb|)=i £ «D(|a)(b|) ® |a)(b|. (9) 

a,be{i,j} a,be{i,j} 

Since O maps pure states to states that are nearly pure, we know that the purity of p satisfies 
tr(p 2 ) ^ (1 — e) 2 ^ 1 — 2e. Evaluating the purity using Equation (O gives 

l-2 £ ^tr(p 2 ) = i(trO(|i)(i|) 2 + trO(|j)(j|) 2 + 2trO(|i)(j|)(D(|j)(l|)) 

<5 + ^tr((<D(|i)(j|)0(|i)(jir). (10) 

Interpreting the expression trXX* as the sum of the squared singular values of X, Equation ((TO)) 
implies that the operator has largest singular value at least 1 — 4e. Since the sum of the 

singular values of this operator cannot exceed one (as the trace norm does not increase under the 
application of a channel), this implies that it can be decomposed as 

<D(|i}0l) = (l-4e)|(b i )((J) j |+4eY / (11) 

where 14k), |(bi) € 0C are pure and Y is a linear operator on % with || Y|| tr = 1. It remains to show 
that the vectors |4k) and |4>j) are, up to a phase, approximately equal to the vectors |i|k) and 
|i|)j) defined in Equation ©. To do this, we consider the action of O on (|i) + y/2. Since O 
is an e-isometry, the output of O on this state is within trace distance 2e of some pure state |y). 
Combining Equations ^ and ((11) and applying the triangle inequality yields 



IyXtI - \ (Nu>MhI + |(b t >((b,| + ItbjXcbtl + I^X^I) 



^ 5c. 

tr 



Since |y) is pure, for some phases Ci and Cj we have |||(bi}({|)j| — c^.c* It^i.) (^P j I lltr ^ 5e, which 
in turn implies that ||0(|i)(j|) — CtC?|ipi)(ipj|||tr ^ 9e, using Equation (jlT) . Finally, since this is 
true for any i + j, and the case of i = j is Equation (|8]), the previous equation implies that 
maxp || O(p) — ApA* || tr ^ 9e, where A is the operator defined by A|i) = Ct for all i. 

It remains only to show that the operator A ® A preserves symmetric states. To see this, take 
+ |ji) an arbitrary basis element of the symmetric subspace of 'K® 1 . By a simple calculation 

(A ® A) (|tj) + |jt) ) = CiCj Mh> ® |i|>j) + c t Cj ) ® |^i) 7 

which remains invariant under swapping the two spaces. By linearity, conjugation by A <2> Igj also 
preserves the symmetry of states on (IK ® "X]® 1 . It follows that 4> preserves symmetry up to an 
error of 9 e in the trace distance. This implies that the swap test on the output of O (g) O applied to 
a symmetric state returns the symmetric outcome with probability at least 1 — 9e. □ □ 

This theorem shows that Non-isometry £ is in QMA for any constant e satisfying (1 — e)/2 > 9e. 
Together with the QMA-hardness of the problem shown in Theorem H] this gives the main result. 

Corollary 8. For any constant e < 1/19, Non-isometry £ is QMA-complete. 

This also implies that problem of computing the channel purity, as defined by Zanardi and Li- 



dar l20h , over the whole input space is QMA-complete. 



10 



6 Conclusion 



We have shown the computational intractability of the problem of detecting when a quantum chan- 
nel is far from an isometry, or equivalently, when a channel can be made to output a highly mixed 
state. These results show that it is extremely difficult to characterize the worst-case behaviour of 
a quantum computation. This is similar to the classical case, where the problem of determining if 
a circuit can produce a specific output is known to be intractable. 

We have also added to the short but growing list of problems that are known to be complete 
for the complexity class QMA. The Non-isometry problem provides a new way to study this class, 
as it exactly characterizes the difficulty of the problems in the class. It is hoped that this will lead 
to new results about the power of this model of computation. 
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